Tailoring diffusion in analogue spacetimes 
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Diffusive transport is characterized by the scaling law (length)^ oc (time). In this paper we show that 
this relationship is significantly altered in curved analogue spacetimes. This circumstance provides 
an opportunity to tailor diflusion: by a suitable design of the analogue metric, it is possible to 
create materials where diffusion is either faster or slower than in normal media, as desired. This 
prediction can in principle be tested experimentally with optical analogues, curved graphene sheets, 
etc. - indeed with any analogue spacetime. 
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I. INTRODUCTION 

Diffusion phenomena are ubiquitous. They underlie 
thermal and electrical conduction, turbidity, osmosis, 
breathing, but also sound propagation in rooms [3,0 
and streets opinion spreading in crowded halls [3| - 
and so forth. In these sundry systems, the characteristic 
feature of the diffusive behaviour is the diffusion scaling 
law: the spatial spreading of a local impulse grows like 
the square root of time. 

Does this scaling still hold in a curved spacetime? At 
first sight, this question may seem rather academic, given 
that spacetime curvature is strong only in extreme as- 
trophysical (dense stars) and cosmological (early uni- 
verse) situations. We believe it is not. As is now well- 
understood, a variety of materials effectively behave as 
analogue spacetimes. This means that, as far as low 
energy excitations are concerned, the effect of micro- 
scopic interactions can be accounted for by a non-trivial 
Lorentzian metric. Examples of such materials include 
corrugated graphene sheets Q, Bose-Einstein conden- 
sates Jl], metamaterials 0, quantum liquids Q, etc. 
(See @ for an updated review of analogue gravity.) To 
study transport phenomena in these systems, it is im- 
portant to have a good grasp on the effect of spacetime 
curvature on diffusion, and notably on the diffusion scal- 
ing law. 

One way to address this problem is via Einstein's 
stochastic approach to diffusion, exposed in his 1905 pa- 
per on Brownian motion [l^. Walking in these footsteps, 
we derived in [ll| the equation governing the probability 
density of Brownian motion in the effective geometry 

ds'^ = -N^it, x)dt^ + qab{t, ^)da;°c^a;^ (1) 

where (t, x) are coordinates comoving with the medium, 
N is the so-called lapse function^ and qab a D- 
dimensional Riemannian metric describing the intrin- 
sic geometry of the spatial slices. We found that, if 
<lab ~ <lab{x) is static, the generalized diffusion equation 



reads [lT| 

dtP^nl\q{Np). (2) 

Here p is the probability density, n the diffusivity and 
Ag the Laplace-Beltrami operator associated to the spa- 
tial metric qab (see below for definitions). Using this 
equation, it is only a computational matter to obtain 
the mean squared displacement (MSD) as a function of 
time, and therefore to obtain the curvature corrections 
to the diffusion scaling law. It is the purpose of this 
letter to present our solution to this problem. 

What makes this computation particularly interesting 
is the fact that the metric coefficients N and qab are of- 
ten tunable in gravitational analogues. For instance, in 
optical media the lapse function N is nothing but the 
inverse refractive index, which recent techniques (meta- 
materials, non-linear Kerr effect) allow to design very ef- 
ficiently. This circumstance suggests that diffusive trans- 
port can perhaps be tailored in analogue spacetimes: ei- 
ther enhanced or slowed down - at will. By analogy 
with the current research on metamaterials [7|, where 
the propagation of electromagnetic radiation can be tai- 
lored with an effective spacetime metric, one can per- 
haps describe our result as establishing the possibility of 
"meifldi/fusioTi" . 

For clarity, we will investigate this issue by considering 
separately the possibility of temporal and spatial tailor- 
ing of the MSD, viz. by distinguishing the cases where 
N and q are only functions of time (resp. space). Once 
these two cases are understood, it is an easy matter to 
treat the general situation, where N and q depend on 
both time and space. 

Among the possible behaviors for the MSD as a func- 
tion of time, we will find: 

• a diffusive-to-ballistic crossover in hyperbolic spa- 
tial geometries, 

• an exponential growth with a parabolic lapse pro- 
file ("Maxwell's fish eye"), 

• a finite limit as i — >■ oo in the presence of an ana- 
logue event horizon. 
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any function f{t) such that /(O) = with a time- 
dependent lapse function. 
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At the very least, this variety of behaviors demonstrates 
that diffusion in curved spacetimes is an interesting topic 
in its own right. We will comment on one possible ap- 
plication in the conclusion. 



II. DEFINITIONS AND ASSUMPTIONS 

Consider a diffusion process within an in irrotational 
fluid flow. In the instantaneous rest frame (Lagrangian, 
or comoving coordinates), the effective spacetime metric 
can be written as 



-N'^it, x)dt^ + qab{x)dx'' dx^ 



(3) 



where N is the lapse function^ and qab a static D- 
dimensional Riemannian metric describing the intrinsic 
geometry of the constant-i spatial slices E^. 

Fix a spatial point xq G Eq, and consider the Green 
function (heat kernel) ift(a;,a;o) of the generalized diffu- 
sion equation ([5]), viz. the solution with initial condition 



Imi Kt{x,XQ) 



(5(x,Xo), 



(4) 



where 8{x,xq) is the Dirac distribution on the t = 
spatial slice Eq. Hereafter, we will denote {T,(p)t the 
pairing between a distribution T and a test function 
on Ej, so that e.g. {S, (l))t = (pit, xq). 

To avoid dealing with drift effects, we will assume fur- 
thermore that Qab is spherically symmetric about xo, viz. 



qab{p)dx''dx'' = a^{p)dp^ + p^dQ 



D-l- 



(5) 



Here the radial coordinate p has the interpretation p = 
{A/AttY^^ , where A is the area of the t = const., p = 
const, surface, and dfljj_^ is the metric on a unit {D — 1)- 
sphere centered on xo- Without loss of generality, we 
take iV(0,a;o) = I. 

With these assumptions, the squared distance between 
x and xq at time t is given by 



d'^{x,Xo) = / dp' a{p'), 

and we can define the MSQ by 

{d^)t = {Kt,d^)t. 



(6) 



(7) 



In a flat spacetime, where = a = 1, the MSQ is well- 
known to be given by {d'^)t = 2KDt: this is the normal 
diffusion scaling law. In the more general metric con- 
sidered in this paper, however, such a linear behavior is 
valid only in the i — >■ asymptotic regime 



{d% 



2KDt. 



(8) 



At later times, corrections depending on A^ and a are to 
be expected, since the Brownian particle then has had 
enough time to explore its neighborhood, and felt the 
effects of its non-trivial geometry. We call tailoring of 
diffusion the possibility of tuning (d^)* as a function of 
t by means of the parameters N and a. 



III. TEMPORAL TAILORING 

Let us consider first the possibility of temporal tailor- 
ing, by which we mean that A^ — ^{t), and a = 1. In 
this case, the equation for the Green function simplifies 
to 



dtKt - NnAKt, 



(9) 



where A is the standard (time- independent) Laplace op- 
erator. Defining 



s{t) = / dt'N{t') 
Jo 



(10) 



and performing the change of variables 1 1-^ s{t) in equa- 
tion dSl), we get 



(11) 



This is nothing but the standard diffusion equation, 
hence 



{d^)t ^ 2KDs{t). 



(12) 



This result implies that, with a suitably designed A^, one 
can tailor the MSD to any desired function f{t) such 
that /(O) = 0. Indeed, differentiating (d'^)t = f{t) with 
respect to t and using (fT^ . we obtain 



N 



f 
2kD' 



(13) 



A lapse function satisfying this condition will ensure that 
the MSD matches the desired function f{t). 



IV. SPATIAL TAILORING 

Even more intriguing is the possibility of spatial tai- 
loring with a static, inhomogcneous spacetime metric: 

ds^ = -N^{p)dt^ + a^{p)dp^ + p2(if7^_i. (14) 

With this spacetime geometry, we cannot in general com- 
pute the MSD in closed form. However, we can easily 
obtain its i — > asymptotic expansion, as follows. 

From now on, denote L the differential operator 
Ag ( A^ • ) ; so ihaX the generalized diffusion equation @ 
reads 

dtKt = nCKt. (15) 
This equation can be solved formally as 

Kt{x,XQ) = e''"^5(x,XQ), (16) 
or, expanding the exponential in powers of its argument, 

Kt{x, xo) = V ^^/:"<5(a;, xo). (17) 

n=0 
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Now, for any test function we have 

{C6,cj,)t^{5,C^)t = {C^){t,xa). (18) 

where = N/S.q is the formal adjoint of C. Using this 
relationship, we find 

(^')*-E^((^'r^?)(-0'-o). (19) 

T1=0 

Explicitly, the order-n coefficient d\ of the Taylor expan- 
sion {d?)t = J2^=o'^n^" therefore 

dl = ^(NA,{--. {NA, d'))){xo,xo). (20) 

n times 

Finally, observe that the asymptotic expansion (fT^ 
can be resummed formally as 

(d')t = (e''*^>2(xo,xo). (21) 

This means that the MSQ itself can be computed from 
a partial differential equation, namely as ut{xo), where 
Ut{x) is the solution to the equation 

dtut - KC^Ut (22) 

with initial condition uo{x) = d^{x,xo). This last equa- 
tion is referred to as the backward Kolmogorov equation 
associated to ([2]) in the mathematical literature. It can 
be useful to compute the MSD numerically. 



It follows that 



i.e. 




The interpretation of this result is as follows: 

• If the lapse profile is convex (e = 1), diffusion 
speeds up in time, with an exponential growth of 
the MSD with time scale t* a^/2KD. With re- 
spect to normal diffusive scaling, such a profile is 
repulsive. 

• If the lapse profile is concave (e = —1), diffusion 
slows down in time, until it eventually stops at 
t ^ t* = /2kD, when the MSD reaches the finite 
limit cr^. With respect to normal diffusive scaling, 
such a profile is attractive. 

From a relativist's perspective, the stopping of diffusion 
on the surface p = a does not come as a surprise: it 
is an instance of the "freezing" phenomenon typical of 
event horizons. We conjecture that the same behavior 
would occur in the vicinity of any infinite redshift sur- 
face, where the lapse function N vanishes. 



V. TWO EXAMPLES OF SPATIAL TAILORING 

To illustrate these findings, let us consider two cases 
of spatial tailoring of direct relevance for applications: 
a quadratic lapse profile, and a spatial geometry with 
constant curvature. 



A. Parabolic lapse profile 

Assume first that the spatial sections are Euclidean 
(a = 1), and the lapse function is parabolic, 

iV(p) = l + e4- (23) 

Here e ~ ±1 indexes the convexity/concavity of the pro- 
file. Such a profile is used e.g. in graded-index optical 
fibers. If e = 1, it is sometimes referred to as Maxwell's 
fish- eye, and has recently attracted a great deal of atten- 
tion for its perfect lensing properties jl2| . 

From a computational viewpoint, this example is par- 
ticularly simple, because both Ap^ and AA^ are constant, 
given by 2D and 2De/a^ respectively. Thus for n > 1 
we have 

[C^Tp' = Ap^(AiV)"-i = . (24) 



B. Pure spatial curvature 

Another interesting special case is when A^ = 1, viz. 
when the only non-trivial curvature components are spa- 
tial. In this case, already considered in the biophysi- 
cal context [i3| , the order- rt Taylor coefficient of the 
MSD as a function of time (PO)) reduces to 

dl^^A^qd\x^,x^). (27) 

The quantities A^c?^(a;o, xp) are well-known geometric 
invariants, tabulated e.g. in [iBl-^ The first two read 

Aqd^{xo,xo) = 2D, (28) 
Ald^{xa,Xf,) = -^i?,(xo), (29) 

where Rq{xo) is the spatial Ricci curvature. Hence, to 
next to leading order, we find that a negative spatial 
curvature will enhance the MSD, while a positive spatial 
curvature will diminish it. The characteristic time of 
the diffusive-to-ballistic crossover in the hyperbolic case 
is t* = D/nRqixa). 



^ See also for fast numerical algorithms. 
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VI. COMMENTS 

Several comments may be useful to clarify the nature 
of the results discussed in the preceding sections. 

First, let us stress that the main equation of this pa- 
per, namely eq. ([2]), is not merely the standard diffusion 
equation written in a non-standard coordinate system. It 
is a different and more general equation, which reduces 
to the standard diffusion only when the spacetime cur- 
vature and background accceleration both vanish. In the 
presence of spacetime curvature, as for instance in solids 
with disclinations or in gradient-index optical media, it 
leads to genuinely different physical predictions. Fur- 
thermore, the coordinate system used to write ([3]) and 
([2]) is not arbitrary: it is the comoving frame in which 
the background medium is locally at rest. 

Second, the generalized diffusion equation ^ is actu- 
ally well-known in general relativity: it is the Eckart heat 
equation [13. (See [m for more details.) Here thanks 
to the connection with Brownian motion established in 
pj| . it is applied to any diffusion phenomenon in the 
presence of an analogue gravitational field, and not just 
to heat dynamics in dissipative relativistic fluids. 

Third, it is worth emphasizing that the analogue grav- 
ity picture of condensed-matter systems underlying this 
work is only an effective description, which hides a great 
deal of microscopic physics. This approximation is simi- 
lar to the one consisting in describing an optical medium 
by a single scalar field, the refractive index. It is valid 
only when the relevant excitations, as the photons in the 
optical case, have low energy and large wavelength and 
therefore do not couple to microscopic inhomogeneities. 
Thus, we have assumed in particular that transport can 
be described as a local diffusion, involving the Laplace- 
Beltrami operator in space Ag (rather for instance than 
as a diffusion in phase space, as in the Kramers kinetic 
equation). As usual, this implies that the relevant mean 
free path is small compared to other length scales in the 
problem. 

Fourth, the Green function Kt of the generalized diffu- 
sion equation ([2]), unlike its flat-spacetime equivalent, is 
in general not Gaussian. This implies in particular that 
Kt is not fully characterized by the associated MSD. 
For this reason, it will be necessary to push further the 
analysis initiated in this paper in order to get a com- 



prehensive picture of "diffusion tailoring". One way to 
gather further information about the effect of analogue 
gravity on diffusion is to compute the time-evolution of 
the entropy of Kt as a function of the effective geometry; 
this task will be undertaken in a future publication. 



VII. CONCLUSION AND OUTLOOK 

In this paper we have computed the MSD for Brown- 
ian motion in a curved, spherically symmetric, analogue 
spacetime. Depending on the convexity of the lapse func- 
tion and the sign of spatial curvature, we have found 
that the MSD can cither grow faster or slower than in 
flat spacetime. In the extreme case where the analogue 
spacetime geometry presents an infinite redshift surface 
(N ~ 0), diffusion actually stops. In principle, any one 
of the many condensed-matter systems behaving as ana- 
logue spacetimes can be used to test these predictions in 
the lab. 

Diffusion is such a common phenomenon in physics, 
chemistry, etc., that it is impossible to list all the possi- 
ble applications of "metadiffusion" . Let us simply men- 
tion one, related to Anderson localization.-^ According 
to the so-called scaling theory, the critical point of the 
localization transition in a finite sample of size L can be 
estimated by the criterion g{L) ~ 1, where g{L) is the 
ratio of the uncertainty of the energy levels due to the 
boundary conditions to the mean energy spacing. The 
former is inversely proportional to the Thouless time : 
the time needed for a diffusive carrier to reach the bound- 
aries of the sample. In normal situations, this time is 
given hy Tl = k/L^. In a sample with curved effective 
geometry, however, we have found that this relationship 
can be modified, and actually tailored. Can this be used 
to control the localization threshold of analogue space- 
times? This question will be addressed in a future work. 
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Note added in proof. The concept of "diffusion tailor- 
ing" was also discussed in the context of heat dynamics 
in fldl], which appeared after the first submission of this 
article to Phys. Rev. E. 
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